Mathematics  Research  Center 
University  of  Wisconsin — Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 


May  1985 


(Received  May  6,  1985) 


OTiC  FILE  COPJf 


Sponsored  by 

U.  S.  Array  Research  Office 

P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Approved  for  public  release 
Distribution  unlimited 


Wxicre 

fcUG  2  0  ''985 


National  Science  Foundation 
Washington,  DC  20550 


093 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 


DEVELOPMENT  OF  SINGULARITIES  IN  NONLINEAR  VISCOELASTICITY 
J.  A.  Nohel  and  M.  Renardy 

Technical  Summary  Report  #2821 

May  1985 

ABSTRACT 

We  discuss  the  motion  of  nonlinear  viscoelastic  materials  with  fading 
memory  in  one  space  dimension.  We  formulate  the  mathematical  problem,  survey 
results  for  global  existence  of  classical  solution  to  the  initial  value 
problem  if  the  data  are  sufficiently  small,  and  discuss  in  detail  the 


development  of  singularities  in  initially  smooth  solutions  for  large  data. 
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DEVELOPMENT  OF  SINGULARITIES  IN  NONLINEAR  VISCOELASTICITY 


J.  A.  Nohel  and  M.  Renardy 
1 .  Introduction  and  Discussion  of  Results 

In  this  paper  we  discuss'^he  motion  of  nonlinear  viscoelastic  materials 
with  fading  memory  in  one  space  dimension.  We  concentration  viscoelastic 
solids  and  briefly  remark  on  similar  results  for  fluids.  After  formulating 
the  mathematical  problems ,  -we  survey  results  for  global  existence  of  classical 
solutions  to  the  initial  value  problem,  provided  the  initial  data  are 
sufficiently  small.  Me-- then  discuss  in  some  detail  the  development  of 
singularities  in  initially  smooth  solutions  for  large  data. _ - 

We  consider  the  longitudinal  motion  of  a  homogeneous  one- dimensional  body 
occupying  an  interval  B  in  a  reference  configuration  and  having  unit 
reference  density.  For  simple  materials,  the  stress  a  at  a  material  point 
x  is  a  nonlinear  functional  of  the  entire  history  of  the  strain  e  -  ux  at 
the  same  point  x  (here  u  denotes  the  displacement).  In  this  paper,  we 
confine  ourselves  to  the  following  model  problem,  which  can  be  motivated  as  a 
natural  generalization  of  Boltzmann's  constitutive  relation  for  linear 
viscoelasticity  [1]  (the  derivation  of  similar  results  in  a  variety  of  other 
models  will  be  discussed  in  a  later  paper) 


o(x,t)  =  4>(e(x,t))  +  J1  a'  (t-T)\|>(e(x,T)  )dt  , 


(x  e  B,  -»  <  t  <  «) 

Here  $  and  are  given  smooth  functions  R  +  R  with 


(1.1) 


$(0)  =  +(0)  =  0,  >  0,  >  0  , 


(1.2) 


and  for  physical  reasons  the  relaxation  function  a  :  [0,®)  R  is  positive. 
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decreasing,  convex,  and  a'  6  L^[0,°°),  where  '  denotes  the  derivative.  The 
conditions  on  a  imply  that  the  stress  relaxes  as  time  increases  and  that 
deformations  which  occurred  in  the  distant  past  have  less  influence  on  the 
present  stress  than  those  which  oc-  curred  more  recently.  Since  only  a' 
occurs  in  the  equation,  we  may  use  the  normalization  a(®)  =  0.  In  the 
rheological  literature  the  relaxation  function  a  is  often  taken  to  be  a 
finite  linear  combination  of  decaying  exponentials  with  positive  coefficients 
obtained  by  a  least  square  fit  to  experimental  data. 

when  (1.1)  is  substituted  into  the  balance  of  linear  momentum,  the 
following  integrodif f erential  equation  for  the  displacement  u  results 

utt  =  <l»(ux)x  +  a'*^(ux)x  +  f  ,  x  e  B,  t  >  0  .  (1.3) 

Here  *  denotes  the  convolution  (a*0)(t)  =  a(t-T )0(t )dx ,  and  f  is  the 
sum  of  an  external  body  force  and  the  history  term  a'(t— r)\|>(u  (x,t  ) )  dt. 

,  _oo  X  X 

An  appropriate  dynamical  problem  is  to  determine  a  smooth  function 
u  :  B  x  (0,0°)  ♦  R  which  satisfies  (1.3)  together  with  appropriate  boundary 
conditions  if  B  is  bounded  and  which  at  t  =  0  satisfies  prescribed  initial 
conditions 

u( x, 0 )  =  Uq ( x ) ,  u^XjO)  =  u ^  (  x ) ,  x  e  B 

for  certain  smooth  functions  Uq  and  u^.  To  avoid  technical  complications, 
we  assume  in  the  following  that  f  =  0.  We  restrict  the  discussion  to  the 
case  when  B  =  R  and  thus  obtain  the  Cauchy  problem 

utt  =  ♦  a’*,J'<ux)x  '  x  e  R#  t  >  0  ,  (1.4) 

u(x, 0 )  =  Uq ( x ) ,  u^(x,0)  =  u i ( x )  ,  X  e  R  .  (1.5) 

When  a'  =  0  and  $  satisfies  (1.2),  the  body  is  purely  elastic.  In 
this  case  it  is  well  known  (see  Lax  [14],  MacCamy  and  Mizel  [15],  Klainerman 
and  Majda  [13]),  that  in  general  the  Cauchy  problem  (1.4),  (1.5)  does  not  have 
globally  smooth  solutions,  no  matter  how  smooth  and  small  the  initial  data  are 
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chosen.  The  initially  smooth  solution  u  develops  singularities  (shock 
waves}  in  finite  time. 

If  a*  %  0  and  a  satisfies  the  sign  conditions  above,  the  fading  memory 
term  in  (1.4)  introduces  a  weak  dissipation  mechanism.  Significant  insight 
into  the  strength  of  this  mechanism  was  gained  by  the  work  of  Coleman  and 
Gurtin  [2] ,  who  studied  the  growth  and  decay  of  acceleration  waves  in 
materials  with  memory.  They  showed  that  the  amplitude  q(t)  of  an 
acceleration  wave  propagating  into  a  homogeneously  strained  medium  at  rest 
satisfies  a  Bernoulli-Riccati  ordinary  differential  equation.  The  coefficient 
of  q  in  this  equation  is  proportional  to  a  second  order  elastic  modulus, 
which  is  given  by  4>"  in  our  model  problem,  and  there  is  a  linear  damping 
term  proportional  to  a'(0).  Thus  amplitude  q(t)  =  [utfc]  decays  to  zero 
as  t  ♦  <">,  provided  |q(0)|  is  sufficiently  small.  On  the  other  hand,  if 
<t>"  ^  0,  then  q(t)  *  *  in  finite  time  if  |q(0)|  is  large  enough,  and 
q(0)  is  of  a  certain  sign. 

This  suggests  that,  under  appropriate  assumptions  on  ^  and  a,  the 
Cauchy  problem  (1.4),  (1.5)  should  have  globally  defined  classical  (C2) 
solutions  for  sufficiently  smooth  and  small  initial  data  uQ,  u-,,  while 
smooth  solutions  should  develop  singularities  in  finite  time  if  the  initial 
data  are  large  in  an  appropriate  sense.  Such  a  global  existence  result  for 
small  data  was  recently  established  by  Hrusa  and  Nohel  [10]  using  delicate  a 
priori  estimates  obtained  by  combining  an  energy  method  with  properties  of 
Volterra  equations  (even  in  the  presence  of  a  small  body  force).  We  refer  to 
a  recent  survey  [9]  for  earlier  small  data  results  on  initial  boundary  value 
problems  modelling  the  motion  of  finite  viscoelastic  bodies,  and  for  technical 
simplif ications  of  the  analysis  in  the  special  cases  $  =  or  a(t)  =  e-t. 

For  the  global  results  the  Cauchy  problem  is  more  difficult  than  the  finite 


body  problem  because  the  Poincare  inequality  is  not  available  to  estimate 


lower  order  derivatives  from  higher  order  derivatives. 

The  remainder  of  our  discussion  will  deal  with  the  formation  of 
singularities  in  finite  time  from  smooth  solutions  of  the  Cauchy  problem 
(1.4),  (1.5).  For  the  special  case  <p  =  ty,  Markowich  and  Renardy  [17]  have 
obtained  numerical  evidence  for  the  formation  of  shock  fronts  in  finite  time 
from  large  data,  and  Hattori  [7]  has  shown  that,  if  <j>"  %  0  and  if  the  body 
B  is  finite,  then  there  are  smooth  initial  data  (which  he  does  not 
characterize)  for  which  the  corresponding  Dirichlet-initial  value  problem  does 
not  have  a  globally  defined  smooth  solution.  On  the  other  hand,  Hrusa  [8]  has 
shown  that  if  $  is  linear  and  only  i|/  is  allowed  to  be  nonlinear,  then  the 
Cauchy  problem  (1.4),  (1.5)  does  have  globally  smooth  solutions,  even  for 
large  smooth  data.  Therefore,  we  shall  restrict  ourselves  to  the  case  when 
?  0,  at  least  over  the  range  of  the  solution.  The  case  when  41"  changes 
sign  will  require  further  refinements. 

An  important  ingredient  in  the  analysis  (which  is  also  important  for  the 
global  theory)  is  the  following  local  existence  result  which  is  established  by 
combining  Banach's  fixed  point  theorem  on  an  appropriate  function  space  with 
standard  energy  estimates  and  Sobolev's  embedding  theorem. 

Proposition  1 ; 

Assume  that  <J>,  t  e  C3(R)  satisfy  (1.2);  assume  a,  a', 

1  ( * ) 

a"  e  L„  [0,®),  and  there  is  a  constant  ic  >  0  such  that 

roc  -  - 

4>'(5)  >  K,  5  e  R  . 


<*) 

Here  the  square  bracket  means  integrability  up  to  0.  No  sign  condition  on  a 
are  required,  but  a'(0)  finite  is  essential. 


2  9 

Assume  that  Uq  e  L^^R)  and  that  u^,  u^  6  H  (R).  Then  the  Cauchy  problem 
(1.4),  (1.5)  has  a  unique  classical  solution  u  e  C2(R  x  [0,TQ) )  defined  on 
a  maximal  interval  (Q,Tq).  _If_  Tq  is  finite,  then 

sup  ^xx**'*5!  +  ^xt**'*^1  ="  ' 

Rx[0,T0) 

The  proof  of  Proposition  1  is  almost  identical  to  that  of  Theorem  2.1  of 
[6],  and  we  omit  the  details;  only  certain  readily  available  energy  estimates 
for  lower  order  derivatives  are  needed.  The  characterization  of  the  maximal 
interval  of  existence  is  established  by  combining  the  energy  estimates 
obtained  in  [6]  with  a  Gronwall  inequality  argument.  We  remark  that  the 
energy  estimates  used  in  the  proof  of  Proposition  1  yield  time -dependent 
bounds  which  cannot  be  used  to  obtain  global  estimates.  These  can  only  be 
constructed  by  taking  advantage  of  the  damping  mechanism  induced  by  the  memory 
term  under  appropriate  sign  conditions  on  a  and  by  assuming  the  initial  data 
to  be  small  (see  [10]  for  details). 

The  assumptions  concerning  the  kernel  a  in  Proposition  1  imply  that 
a'  is  absolutely  continuous  on  [O,00).  Recently,  Hrusa  and  Renardy  [11] 
established  a  result  similar  to  proposition  1  (and  proved  a  global  existence 
result  for  small  data  for  bounded  bodies)  under  assumptions  which  permit  a 
singularity  in  a'  at  t  =  0  (e.g.  a'(t)  -  -t“  \  0  <  a  <  1  as  t  0+). 
Such  singularities  are  relevant  for  certain  popular  models  of  viscoelastic 
materials. 

Our  main  result  on  development  of  singularities  for  large  enough  data  is 

theorem  1: 

o  1 

Let  $,  e  CJ(R)  satisfy  (1.2)  and  assume  a,  a',  a"  e  L£ocf°'“)* 

Assume  that  $"(0)  /  0.  Then,  for  every  T^  >  0,  we  can  choose  initial 
data  ufo,  u1  e  C^(R)n  L°°(R)  such  that  the  maximal  time  interval  of 


-5- 


existence,  given  by  Proposition  1,  for  the  smooth  solution  of  the  Cauchy 

problem  (1.4),  (1.5)  cannot  exceed  T, .  More  precisely,  if  sup  |u„(x)|  and 

xeR 

sup  [ u  ( x )  |  are  sufficiently  small,  while  u'q  ( x )  and  u  ^  ( x )  assume 
x6R  ' 

sufficiently  large  values  with  appropriate  signs,  then  there  is  some 
t  <  T1  such  that 

sup  #  {|uxx(x,t)|  +  |uxt(x,t)|}  =  »  ,  (1.6) 

Rx[0,t  ) 

while 

sup  *  {|u  (x,t)|  +  |u  ( x, t ) | }  <  «  (1.7) 

Rx[0,t  ) 

(and  in  fact,  this  latter  quantity  remains  small). 

In  view  of  the  analogy  with  hyperbolic  conservation  laws  and  the 

numerical  evidence  [17],  it  is  to  be  expected  that  a  blow-up  as  established  by 

Theorem  1  will  lead  to  the  development  of  a  shock  front. 

The  method  of  the  proof,  sketched  in  section  2,  is  to  show  that  the 

memory  term  is  in  fact  of  lower  order  than  the  elastic  term  <MUX^X  and  can 

be  treated  as  a  perturbation.  While  considerably  more  technical,  the  proof  is 

a  generalization  of  the  approach  of  Lax  [14]  for  showing  the  development  of 

singularities  for  the  quasilinear  wave  equation 

u,.,.  =  d>(u  ) 
tt  x'x 

Theorem  1  was  established  independently  by  Dafermos  [4]  using  an  approach 
which  is  different  from  ours  but  similar  in  spirit.  The  result  can  also  be 
established  by  modifying  the  results  of  F.  John  [12]  and  extending  them  to 
systems  of  quasilinear  hyperbolic  conservation  laws  which  contain  lower  order 
source  terms  (F.  John,  private  communications). 

Similar  results  for  first  order  model  problems  were  derived  by 
Malek-Mandani  and  Nohel  [16]  and,  using  different  methods,  by  Renardy  [18]  and 


Dafermos  [3] 


A  particular  case  of  the  model  equation  studied  in  this  paper  leads  to  a 


model  fOi.  shearing  flows  of  viscoelastic  fluids  studied  recently  by  Slemrod 
[ 20] -  With  v(x,t)  denoting  the  velocity  of  the  fluid  in  simple  shear, 
Slemrod  studies  the  problem 

vt  =  a*<Mvx)x  >  (x  e  R,  t  >  0)  , 


v( x, 0 )  =  vQ(x)  ,  (x  €  R) 


(1-8) 


for  the  special  case  a  =  e-t.  Problem  (1.8)  leads  to  a  Cauchy  problem  of  the 
form  (1.4),  (1.5)  after  differentiation  with  respect  to  time.  Thus  Theorem  1 
can  be  used  to  get  a  blow-up  result  for  this  problem,  like  the  result  found  by 
Slemrod  for  a(t)  =  e-t.  The  global  existence  of  solutions  for  small  data 
follows  from  [5,  Theorem  4.1].  Other  popular  models  for  viscoelastic  fluids 
have  been  analyzed  by  the  method  used  in  this  paper;  the  results  will  be 
published  elsewhere. 


2.  Develops*  f  t  of  Shocks 

In  this  section,  we  sketch  the  proof  of  Theorem  1  establishing  the 
development  of  shocks  from  initially  smooth  solutions  of  the  Cauchy  problem 
(1.4),  (1.5)  in  finite  time.  For  simplicity,  most  of  the  analysis  will  be 
carried  out  for  the  special  case  a(t)  =  e_t-  the  proof  for  more  general 
relaxation  functions  as  well  as  for  a  more  general  class  of  model  equations 
will  be  carried  out  in  a  forthcoming  paper. 

We  begin  by  transforming  (1.4)  to  an  equivalent  system.  We  let 
w  =  ux,  v  =  ut,  and  write  the  constitutive  assumption  (1.1)  in  the  form 

o  =  <)>(w)  -  z  ,  z  =  -a'*ty(w) 

Since  we  have  assumed  >0,  the  first  of  these  equations  can  be 

solved  for  w, 

w  =  4>  ^(o+z)  =:  g(a.z)  , 
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and  g  is  a  smooth  function  of  o  e  R,  z  6  R.  As  long  as  the  solution 
remains  smooth,  the  Cauchy  problem  (1.4),  (1.5)  is  equivalent  to  the  first 
order  system 


ot  =  c2(o,z)v^  +  a*  (0)\J>(g(cs,z) )  +  a"*\|)(  g  ( a ,  z ) )  , 
zt  =  -a'  (O)i|*(g(0,z) )  -  a"*i|;(g(c,z) ) 


The  initial  conditions  become 

v(x,0)  =  u-,(x),  o(x,0)  =  <t>(Ug(x)),  z ( x , 0 )  =  0  .  (2.2) 

By  c  we  have  denoted  the  wave  speed 

1/, 

c  (a,  z)  :=  [  <t>  *  (g(0,z)  )]'2  ; 

c  is  a  smooth  function  of  o  and  z.  The  system  (2.1)  is  hyperbolic,  and 
its  eigenvalues  are  +c,  -c  and  0.  Under  the  assumptions  of  Proposition  1, 
a  ^-solution  exists  on  some  maximal  interval  R  x  [0,TQ).  if  Tq  is  finite, 
then  v,  a,  z  or  one  ot  their  first  derivatives  must  become  infinite  as 
t  *  Tn.  It  is  immediate  from  equation  (2.1)  that  o  ,  a  ,  z*.  and  zv  will 
remain  bounded  as  long  as  v,  o,  z,  vfc  and  vx  are  bounded. 

To  proceed  further,  we  extend  the  classical  approach  of  Lax  [14]  for 
first  order  hyperbolic  2  x  2-systems.  We  define  "approximate”  Riemann  in¬ 
variants  by  those  quantities  which  would  be  the  classical  Riemann  invariants 
if  z  in  the  first  two  equations  of  (2.1)  were  treated  as  a  parameter.  These 
quantities  are  given  by 


r  =  r  ( v ,  o,  z) 
s  =  s(v,o,z) 

$(a,z)  = 

0 


=  v  +  $(o, z) 
=  v  -  $>(o,z) 

. 

c(C, z) 


(2.3) 
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0. 


thout  loss  of  generality  we  may  take  o = 


Since  $>^(o,z) 


1 


a,z ) 


>  0,  this  correspondence  is  smoothly  invertible,  and  we  have 


r+s  %  r-s 

v  =  ——  ,  4>(o,z)  =  — 


i  the  following,  we  assume  a(t)  =  e_t-  Then  (2.1)  takes  the  simple  form 


v.  =  a 
t  x 


o  =  c  ( a , z ) v^  -  i|»(  g  (  o  ,  z )  )  +  z  , 


(2.4) 


zfc  =  (Hg(o,z)  )  -  z 


e  now  differentiate  r  and  s  along  the  c  and  -c  characteristics 
espectiv^ly,  and  z  along  the  zero  characteristic  (i.e.  we  form 
t  -  crx,  st  +  csx  and  zt).  This  leads  to  the  following  first  order 
lyperbolic  system  equivalent  to  (2.4),  (2.2) 

r.  -  Ar  =  -Bz  +  CD  , 
t  x  x 

S^  +  As  =  -Bz  -  CD  , 

t  X  X 


(2.5) 


z 


t 


=  0  , 


fith  the  initial  data 


r  (  x  ,  0  )  =  u  ^  ( x  )  +  4>(<(i(Uq(x)  ),D)  , 

s  (  x  ,  0 )  =  u  ^  ( x  )  -  $(<j>(u^(x)  )  ,0)  ,  (2.6) 

z  (  x ,  0  )  =  0  ; 


-y- 


=  A(r,s,z)  :=  c ( a ( r , s , z ) , z )  >  0  , 


=  B(  r ,  s ,  z  ) 


:=  c( a(r , s, z ) , z)$2 ( o(r , s, z) , z )  , 


=  C(r,s,z) 


<5  (o(r,s,z),z)  -  — — - r - r  , 

z  c(o(r,s,z),z) 


=  D(  r  ,  s  ,  z  )  :=  \J>(g(0(r,s,z)  ,z)  )  -  z 


To  establish  the  development  of  shocks  in  finite  time,  we  study  the 
lution  along  characteristics  of  the  quantities 


v  +  -  , 

x  c(o,z) 


v  -  -  , 

X  c( 0, z) 


(2.8) 


zx-  Note  that  if  z  were  a  constant  parameter,  then  p  and  t  would 
the  x-deri vati ves  of  r  and  s.  We  have  vx  =  (  p+T ) ,  ax  ~  C*p-T)'  and 


/  2 

[c  )  (0, z)  =  2cc 


=  4>"(g(o*z> ) 

0  4> ’  (g(  0,  Z) ) ' 


edious  but  straightforward  calculation  using  the  relations  (obtained  by 
forentiating  (2.4)) 


Vtx  °xx  ' 

o  =  c2(o,z)v  +  ;c2]  (o,z)a  v 

tx  XX  '  '0  XX 

+lc2](o,z)zv  -D  , 

Z  XX  X 


(2.9) 


ztx  Dx  ' 


1 1  is  the  system 


kt  “  C>'x  =  - 4 -  p(p-T)  +  °UpMzxI  +  I  T  I  I  zx  i 

+  |p I  +  M  +  |zxl)  , 


t  x 


- -  T(P-T)  +  0(  [  P  |  jzj  +  I  T  |  |  -x  | 


(2.10) 


+  p  +  T  +  Z. 


u(  lei  V  I  r I  +  I z  , 


-In- 


subject  to  the  initial  data 


p(x,0)  =  u'^tx)  +  4>*  (u^(x)  )^2  ug(x)  , 

t(x,0)  =  u'(x)  -  <fr,(ug(x))1/'2u£(x)  , 
zx(x,0)  =  0  . 


(2.11) 


The  cross  product  terms  pt  in  (2.10)  are  eliminated  if  one  considers  the 

Vo  Vo 

characteristic  derivatives  of  c (o,z)  *p  and  c(o,z)  (see  Lax  [14]  and 

Slemrod  [19]).  We  find 

“  c  |j){c^2p>  =  Y(c/2p)2  +  0(  |p|  |zj  +  |t|  |zj 

+  |p|  +  | T |  +  |zx|)  , 

(lt  +  C  4t,(c1/2t)  =  V(c/2t)2  +  °( ( P |  | *x  I  +  Mlzxl  (2.12) 

+|p|+lTl+lzxl>  » 
zxt  =  0(lpl  +  M  +  lzxl )  • 

Here  the  coefficient  function  y  is  given  by 


Y  -  Y<«.»)  ■  1  . 

For  definiteness,  let  us  assume  $"(0)  >  0  (the  discussion  for  <j>"(0)  <  0  is 
analogous).  We  take  initial  data  with  the  following  properties:  ug  and 
u^  (and  hence  p(x,0),  t(x,0)  as  well  as  z(x,0)  =  0)  are  uniformly  small, 
and  p(x,0),  t(x,0)  are  such  that  at  least  one  of  them  has  a  large  positive 
maximum  (by  choosing  ug  or  u}  or  both  sufficiently  large).  At  the  same 
time,  the  maxima  of  -p  and  -t  should  not  be  too  large. 

As  long  as  (r,s,z)  remains  within  a  given  neighborhood  U  of  0,  we 
have  upper  and  lower  bounds  for  the  coefficients  occuring  in  (2.12),  in 
particular,  we  have  a  positive  lower  bound  y_  for  y.  We  shall  see  later 


that  (r,s,z)  will  remain  in  U  up  to  the  time  of  blow-up  if  they  are  small 


enough  initially  and  if  we  make  the  maximum  of  p(x,0)  or  t(x,0)  large 
enough . 

For  every  t  >  0,  we  now  set 

h(t)  =  maxfmax  p (..,£),  max  x(x,t)] 
x  x 

From  (2.12),  we  find  that,  as  long  as  (r,s,z)  e  U,  while  h(t)  is  large 

and  maxjz  |  «  h(t),  we  have,  for  some  positive  constants  Yq  and  tc 

x 

(lt>+  h(t)  >^0(h(t))2'  and  max!zxtl  *  ,ch(t)  <<  (h(t))2  • 

X 

Initially,  we  have  |zx|  =  0  and  it  follows  from  these  inequalities  that  it 

will  remain  small  compared  to  h(t).  We  also  find  that  h(t)  becomes 

infinite  in  finite  time.  Since  there  is  also  some  constant  such  that 

(■§--)  h(t)  <  Y..(h(t))#  it  can  be  shown  that,  with  t  denoting  the  blow-up 
at  +  1 

°1  C2 

time  of  h,  we  have  — - —  <  h(t)  <  — - —  for  some  constants  c1  and  c2. 

t  -t  t  -t 

The  third  equation  of  (2.12)  then  implies  that  |zx|  grows  at  most 
logarithmically  as  t  -*■  t*.  Since  log(t*-t)  is  integrable,  equations  (2.5) 
imply  that  r,  s,  and  z  remain  bounded  and  in  fact  small  if  their  initial 
data  are  small,  and  t*  is  small  (which  is  the  case  if  h(0)  is  large).  In 
this  way,  we  can  choose  the  data  such  that  (r,s,z)  will  in  fact  remain  in 
U  up  to  the  time  of  blow-up.  This  completes  the  sketch  of  the  proof. 
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